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THE  VARIATIONAL  METHOD  FOR  AERODYNAMIC  OPTIMIZATION 
USING  THE  NAVIER-STOKES  EQUATIONS* 

BAMBANG  L  SOBMARWOTQt 

Abstract.  This  report  describes  the  formulation  of  an  aerodynamic  shape  design  methodology  using  a 
compressible  viscous  flow  model  based  on  the  Reynolds- Averaged  Navier-Stokes  equations.  The  aerodynamic 
shape  is  described  by  a  set  of  geometrical  design  variables.  The  design  problem  is  formulated  as  an  optimiza¬ 
tion  problem  stated  in  terms  of  an  aerodynamic  objective  functional  which  has  to  be  minimized.  The  design 
scheme  employs  a  gradient-based  optimization  algorithm  in  order  to  obtain  the  optimum  values  of  the  design 
variables.  The  gradient  of  the  aerodynamic  functional  with  respect  to  the  design  variables  is  computed  by 
means  of  the  variational  method,  which  requires  the  solution  of  an  adjoint  problem.  The  formulation  of  the 
adjoint  problem  is  described  which  leads  to  a  set  of  adjoint  equations  and  boundary  conditions.  Using  the 
flow  variables  and  the  adjoint  variables,  an  expression  for  the  gradient  has  been  constructed.  Computational 
results  are  presented  for  an  inverse  problem  of  an  airfoil.  It  will  be  shown  that,  starting  from  an  initial 
geometry  of  the  NACA  0012  airfoil,  the  target  pressure  distribution  and  geometry  of  a  best-fit  of  the  RAE 
2822  airfoil  in  a  transonic  flow  condition  has  been  reconstructed  successfully. 

Key  words,  aerodynamic  optimization,  airfoil  design,  variational  method,  optimal  control,  inverse 
design,  Navier-Stokes  equations 

Subject  classification.  Applied  and  Numerical  Mathematics 

1.  Introduction.  Methodologies  for  solving  aerodynamic  shape  design  problems  can  be  distinguished 
into  two  classes:  (i)  inverse  methodology  and  (ii)  optimization  methodology.  The  distinction  is  based  on 
how  the  design  problem  is  formulated. 

In  the  inverse  methodology,  the  design  problem  is  posed  in  terms  of  a  prescribed  target  pressure  distri¬ 
bution  which  has  to  be  reahzed  on  the  surface  of  the  shape.  The  designer  is  assumed  to  be  able  to  prescribe 
the  target  pressure  distribution  in  such  a  way  that  it  reflects  required  aerodynamic  characteristics  like  lift, 
drag,  pitching  moment,  and  boundary  layer  properties  which  determine  the  aerodynamic  performance.  In¬ 
verse  methods  assist  the  designer  by  constructing  an  aerodynamic  shape  which  generates  the  target  pressure 
distribution  (Refe.  [22],  [11],  [10],  [5]). 

In  the  optimization  methodology,  the  design  problem  is  posed  as  a  minimization  problem  of  an  aerody¬ 
namic  objective  functional  subject  to  constraints  on  the  geometric  and  aerodynamic  properties.  Optimization 
methods  assist  the  designer  in  locating  the  minimum  of  the  objective  while  satisfying  the  constraints.  Prom 
the  practical  point  of  view,  aerodynamic  optimization  methods,  pioneered  by  Hicks  et  al.  [14],  are  more 
attractive  since  these  methods  can  handle  a  large  class  of  design  problems,  including  those  classified  as  in¬ 
verse  problems.  This  report  describes  a  contribution  to  the  development  in  the  aerodynamic  optimization 
methodology. 

Aerodynamic  optimization  methods  can  be  distinguished  into  two  categories:  (i)  global  methods  and  (ii) 
local  methods.  Global  methods,  such  as  those  based  on  the  geneMc  algorithm  [9],  are  aimed  at  obtaining  the 
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global  optimum.  These  methods  are  most  useful  for  cases  in  which  multiple  minima  are  present  in  the  design 
space.  It  is  widely  known,  however,  that  global  methods  incur  large  computational  effort,  where  hundreds 
or  even  thousands  of  flow  analyses  may  be  needed  before  the  global  optimum  can  be  found. 

Local  methods  use  the  information  on  the  gradient  of  the  objective  for  locating  the  optimum.  Therefore, 
for  cases  with  multiple  minima,  local  methods  are  limited  to  produce  only  one  of  the  minima  (i.e.,  the  local 
optimum),  the  actual  value  of  which  depends  on  the  starting  point  of  the  optimization  process.  Because 
of  the  modest  computational  requirement,  and  since  any  local  optimum  represents  an  improvement  over  an 
existing  design,  local  methods  are  very  useful  design  tools.  The  method  described  in  this  report  belongs  to 
this  category. 

Recent  developments  in  the  gradient-based  aerodynamic  optimization  methodology  suggest  that  two 
tnam  streams  may  be  distinguished;  (i)  the  method  of  sensitivity  analysis  and  (ii)  the  variational  method. 
This  distinction  is  based  on  how  the  gradient  is  computed. 

The  formulation  of  the  gradient  using  the  sensitivity  analysis  method  (Refs.  [27],  [24],  [15]),  is  done  on 
a  discrete  level,  which  means  that  one  must  deal  with  the  discrete  form  of  the  flow  equations.  This  method 
has  the  advantage  that  the  sensitivities  of  the  flow  properties  on  the  grid  points  can  be  determined.  Once 
these  become  available,  the  gradient  of  an  aerodynamic  fimctional  can  be  computed  easily  using  the  chain 
rule.  However,  the  computational  effort  strongly  depends  on  the  number  of  design  variables.  For  each  design 
variable,  a  sensitivity  equation  in  the  form  of  a  (large)  linear  system  of  equations  must  be  solved.  The 
computational  requirement  can  therefore  be  prohibitive  if  a  large  design  space  is  to  be  covered. 

The  formulation  of  the  gradient  using  the  variational  method  can  be  done  either  on  a  discrete  level 
(Refs.  [2],  [23],  [6],  [21])  or  a  continous  level  (Refs.  [25],  [18],  [20],  [19],  [26],  [1],  [17],  [16],  [28]).  This  method 
needs  the  values  of  the  so-called  adjoint  variables  as  the  solution  of  a  set  of  adjoint  equations.  The  numerical 
solution  procedure  for  solving  the  flow  equations  can  be  adopted  for  solving  the  adjoint  equations.  The 
gradient  is  expressed  in  terms  of  the  flow  variables  and  the  adjoint  variables.  The  computational  effort  for 
obtaining  the  gradient  is  not  determined  by  the  mimber  of  design  variables.  Instead,  it  is  determined  by  the 
number  of  adjoint  equations  that  must  be  solved,  which  is  equal  to  the  number  of  aerodynamic  functionals  of 
the  aerodynamic  objective  and  constraints.  Anticipating  that  the  number  of  design  variables  is  significantly 
larger  than  the  number  of  aerodynamic  functionals,  which  is  tnie  in  many  practical  cases,  the  variational 
method  has  a  significant  advantage  over  the  method  of  sensitivity  analysis. 

This  report  describes  the  design  approach  utilizing  the  variational  method  for  airfoil  design  using  the 
compressible  Reynolds-Averaged  Navier-Stokes  (BANS)  equations.  Recently,  the  author  [28]  has  demon¬ 
strated  the  feasibility  of  the  approach  in  dealing  with  inverse  problems  and  constrained  drag-reduction 
problems,  where  the  compressible  viscous  flow  model  based  on  the  RANS  equations  was  used  for  the  flow 
calculations.  An  analytical  expression  of  the  adjoint  equations  was  formulated  based  on  the  continous  form 
of  the  aerodynamic  functional  and  the  RANS  equations.  In  the  previous  work,  however,  considerations  from 
the  physics  of  the  boundary  layer  must  still  be  taken  for  obtaining  an  approximation  of  the  gradient,  despite 
the  success  in  obtaining  true  viscous  adjoint  solutions.  Although  the  approximation  can  lead  to  useful  results, 
as  shown  in  Ref.  [28],  it  is  desirable  to  have  a  gradient  expression  which  is  derived  consistently  using  the 
RANS  equations.  The  objective  of  the  present  study  is  therefore  to  obtain  the  true  RANS-based  gradient 
expression. 

2.  Statement  of  the  Design  Problem.  The  design  problem  being  addressed  is  formulated  as  a 
minimization  problem  of  an  aerodynamic  functional 

(2.1)  Minimize  .F(Q,  0) 
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where  Q  is  the  vector  of  flow  variables,  and  0  is  a  vector  representing  the  geometric  parameters  that  define 
the  aerodynamic  shape.  The  vector  $  is  treated  as  the  design  variables,  the  optimal  value  of  which  is  to  be 
determined.  There  is  an  implicit  dependency  of  Q  upon  0  through  the  RANS  equations  for  a  given  onset 
flow  condition. 

Problem  (2.1)  is  to  be  solved  by  means  of  an  iterative  gradient-based  optimization  algorithm.  This 
requires  information  on  the  gradient  of  7  with  respect  to  0  for  each  iterate.  An  efficient  way  of  treating 
the  implicit  dependency  of  Q  upon  0  in  evaluating  the  gradient  is  through  the  variational  method.  In  the 
variational  method,  an  adjoint  problem  must  be  formulated  in  which  a  set  of  adjoint  equations  are  to  be 
solved  subject  to  proper  adjoint  boundary  conditions.  The  gradient  is  expressed  in  terms  of  the  flow  variables 
(i.e.,  the  solution  of  the  flow  problem)  and  the  adjoint  variables  (i.e.,  the  solution  of  the  adjoint  problem). 

3.  The  Reynolds- Averaged  Navier-Stokes  Equations.  Assuming  adiabatic  flow  and  no  external 
forces,  the  time-dependent  RANS  equations  in  the  conservative  form  are  written  as 

(3.1)  ^-HV-F  =  0infl, 

where  fl  is  the  flow  domain,  and  Q  is  the  vector  of  conservative  time-averaged  flow  variables. 


(3.2) 


^  P  ^ 

pu 

pv 

\pEj 


which  are  non-dimensionalized  with  respect  to  the  free  stream.  At  the  steady-state,  equation  (3.1)  becomes 


(3.3)  ,  V  •  F  =  0. 

The  flux  P  consists  of  the  convective,  Fo  and  viscous,  F^,,  flux  vectors, 

(3.4)  F  =  Fc-F^. 


The  convective  flux  vector  is  defined  as 
(3.5) 


Fc  = 


where  and  are  the  Cartesian  components  given  by 

(3.6)  /, 


where  p,  u,  v,  p  and  E  are  the  air  density,  x-  and  y-velocity  components,  pressure,  and  total  energy, 
respectively.  The  viscous  flux  vector  is  defined  as 

fv' 


/  pu  \ 

(  pv  \ 

puu  +  p 

puv 

puv 

>  9c 
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(3.7) 


F.= 


9v 


where  and  are  the  Cartesian  components  given  by 

/  0  \ 

'^xx 


(3.8) 
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Assuming  that  air  behaves  like  a  Newtonian  fluid,  the  elements  Txxj  Tyy  of  the  viscous  stress  tensor 

are  expressed  as 


-4  — ♦  ou 

(3.9)  r**  =  KV  •  V)  +  2/i^, 

-*  -+  dv 

(3.10)  Tyy  =  /(V  •  V)  +  2lX  —  , 


(3.11) 


'^xy  M 


/  dv  du\ 
dy) 


where  I  is  given  by  the  Stokes  hypothesis, 


The  viscosity  y  consists  of  the  dynamic  viscosity  pd  and  the  eddy  viscosity  yt, 


IX  =  /id  +pt, 

where  /td  is  given  in  terms  of  the  onset  flow  condition  by  the  Sutherland’s  law, 

/  T  Too +  110 
/icc  yr^J  T+110’ 

with  T  the  absolute  temperature,  while  /xt  is  defined  by  a  turbulence  model  which,  in  the  present  study,  is 
based  on  the  Baldwin  and  Lomax  model 

The  Cartesian  components  of  the  heat  flux  vector  q  are  defined  by 


dT 

(3.12) 

II 

! 

dT 

(3.13) 
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where  the  thermal  conductivity  coefficient  k  consists  of  the  laminar  part,  and  turbulent  part,  Kt.  These 
are  related  with  the  viscosities  through  the  Prandtl  numbers, 

■n  Md 

Prd-Cp^d> 

Prt  = 

Kt 

with  Cp  the  specific  heat  at  constant  pressure  and  h  the  mass  specific  enthalpy.  The  Prandtl  numbers  are 
assumed  to  have  constant  values  throughout  the  flow,  Prd  =  0.72  and  Pr*  =  0.9,  respectively.  The  total 
energy  E  per  unit  mass  is  defined  as 

E^e+^{u^  +  v^), 

where  e  is  the  internal  energy  per  unit  mass.  The  RANS  equations  are  closed  by  the  equation  of  state  of  a 
calorically  perfect  gas,  given  as 


(3.14) 

P={l-  ^){pE  -  i/>(u2  +  v^)), 

(3.15) 

r=l(£;-lu2  +  n2)). 
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where  '^  =  Cplc^,  with  Cu  the  specific  heat  at  constant  volume.  The  heat  fiuxes  can  be  written  in  terms  of 
the  internal  energy  as 


(3.16) 

(3.17) 


qx  = 

qv  = 


/i  de 

fji  de 
“'’'Pr^’ 


where 

^Pr  ^\Fvd  Pit  J 

On  the  airfoil  surface,  Sa,  the  no-slip  and  adiabatic  boundary  conditions  are  applied.  The  no-slip  boundary 
condition  can  be  expressed  as 

(3.18)  V-n  =  0, 

(3.19)  P-s  =  0, 


where  V  denotes  the  velocity  vector,  while  n  and  s  are  the  unit  normal  and  tangential  vectors,  respectively. 
The  adiabatic  wall  botmdary  condition  reads 

(3.20)  VT  •  n  =  0. 


This  is  formulated  in  terms  of  the  internal  energy  as 

(3.21)  Ve-n  =  0. 


The  boundary  conditions  are  collected  into  a  vector  B  as  follows. 


(3.22) 


^  V  n\ 

B  = 

V-s 

— 

0 

yVe-fiJ 

\o) 

4.  Formulation  of  the  Adjoint  and  Gradient  Equations.  It  is  assumed  that  the  aerodynamic 
functional  :F  takes  the  form  of  a  surface  integral  over  the  airfoil  surface  Sa- 


(4.1) 


dS, 


where  'tp  is  an  exphcit  function  of  the  pressure  p,  the  wall  shear  stress  Tw  and  the  design  variables  0,  with 


Tw  =  p 


d{V-s) 

dn 


The  fimctional  (4.1)  represents  a  large  class  of  design  problems,  including  those  expressed  in  terms  of  lift, 
drag,  and  pitching  moment. 

As  Q  is  obtained  from  the  steady-state  RANS  eqxiations  with  the  boundary  conditions  (3.22),  the 
fimctional  is  independent  of  the  transient  state.  Therefore,  it  is  sufficient  to  consider  the  steady-state 
RANS  equations  (3.3)  and  the  boundary  conditions  (3.22)  in  the  definition  of  a  Lagrangian  C  as  foUows, 


(4.2) 


C  = 


r-BdS, 
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where  A  and  T  are  the  vectors  of  Lagrange  multipliers.  The  Lagrange  multipliers  A,  also  referred  to  as  the 
adjoint  variables,  are  defined  in  Cl  and  consists  of  four  components.  The  Lagrange  multipliers  Y  is  a  vector 
with  three  components  defined  on  Sa¬ 
in  order  to  derive  the  adjoint  and  gradient  equations,  one  must  evaluate  the  variation  of  C,  denoted  as 
6£,  implied  by  the  independent  variations  of  A,  X,  Q,  and  0, 

5C  =  5C\  +  SCt  +  5Cq  +  5Cb- 

The  notation  5Cx  refers  to  the  variation  of  5C  due  to  the  variation  of  A  while  the  other  variables  are  kept 
fixed,  and  similarly  for  5Cr,  etc.  The  variations  5Cx,  6Cr,  and  SCq  are  evaluated  with  0  kept  fixed.  Keeping 
0  fixed  implies  a  fixed  domain  Cl.  For  the  variation  of  A,  X,  and  Q  with  a  fixed  domain  a  prime  notation  is 
introduced  as  A',  X',  and  Q',  respectively. 

4.1.  The  Adjoint  Equation.  The  variation  A'  contributes  to  SC  with 

(4.3)  SCx=  [  A'  (V-F)dn, 

Jn 

which  is  cancelled  by  the  RANS  equations  (3.3).  The  variation  X'  contributes  with 

(4.4)  SCr=  f  r'  BdS, 

J  Sa 

which  vanishes  because  of  the  boundary  conditions  (3.22). 

As  the  RANS  equations  (3.3)  and  the  boundary  conditions  (3.22)  are  satisfied,  giving  SCx  =  0  and 
SCy  =  0,  the  variation  of  C  becomes 

SC  =  SCq  +  SCff. 

The  adjoint  equations  and  boundary  conditions  follow  from  the  condition  that  the  contribution  from  the 
variation  Q'  vanishes,  i.e. 


(4.5)  SCq  =  0. 

The  domain  integral  in  equation  (4.2)  can  be  integrated  by  parts  to  give 

(4.6)  C=  [  i)dS-  f  X-{F-n)dS-f  \-{F-n)dS 

J  Sa  Sa  Szxi 

-  /  F  •  VA  dfl  +  /  r-BdS. 

J  Cl  J  S  a 

The  variation  6Cq  can  be  expressed  as 

SJ^Q=  f  X-iF'-n)dS-f  X-{F'-fi)dS 

-  [  F'  ■yXdCl+  [  r-B'dS, 

Jq  JSa 

where  the  notations  F'  and  B'  refer  to  the  variations  due  to  Q'.  The  flux  vector  F'  can  be  split  into  the 
inviscid  and  the  viscous  part; 

F'  =  F'  -  F' 
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It  is  convenient  to  introduce  the  inviscid  and  viscous  variations,  51  and  5J,  implied  by  and  F^,  respec¬ 
tively,  and  the  variation  5K  implied  by  B',  such  that 


5CQ  =  5rQ  +  51-5J  +  5K. 


The  variation  5Lq  will  be  obtained  with  the  assumption  that 

•  The  variation  of  the  viscosity,  n',  can  be  neglected. 

•  The  variation  of  the  viscous  terms  F:^  on  the  far-field  boundary  Soo  can  be  dropped. 
The  aerodynamic  functional  contributes  with 


/  I  /  1  ,,q 

UTw 


where 


d{V^  •  ^ 


The  inviscid  term  can  be  obtained  as 


6X  =  —  f  (A  •  n)(7  —  l){pEy  dS  —  f  (xi  H - •  n)  dS 

J  Sa  JSa  \  ' 

-  [  {C'^X)  Q'dS-  [{A^  ■  VA)  •  Q' dQ, 

-/fi 


where  A  is  the  Jacobian  of  the  flux  vector  Fc  with  respect  to  Q, 


C  is  the  Jacobian  of  the  normal  flux  defined  on  the  boundaries 

A  is  an  adjoint  velocity  vector  with  the  Cartesian  components  A2  and  A3: 


The  procedure  for  obtaining  dj  is  described  in  the  appendix,  with  the  result  given  as  equation  (A.30).  The 
variation  6IC  can  be  obtained  as 


=  j  TiiY' ■ri)  +  r2{v' ■i}  +  r3^ 


<l_  V{pE)'-n  Vp'-rt 
7  P  (7  - 1)^ 


where  a  is  the  speed  of  sound. 


„  /  P 

“=  a/7-- 
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Substituting  equation  (4.8),  (4.9),  (4.10)  and  (A.30)  into  (4.7)  leads  to 

(4.11)  SCq  =  J^  [^(7 -!)(/>£?)'  + ^4- (A  •«)(7-l)(/»^;y-(Ai  +  :^A4)Ky'.n) 
+(A  •  n)r4  +  (A  • 


-  '(V-A)  +  2;.®5^ 


(F'-n)- 

V{pEy  •  n 

Vp'  •  n  \ 

V 

(7  -  ) 

{V'-n)-p 

d{Xs)  d{X-n) 

dn  ds 

(7  - 1)^/ 

-ff(A.s)  (f'-^ 


+ri(t/'-n)  +  r2(7'-s)  +  r3— 


WjpEY  •  n  V/)^-n 


(7  -  1)P , 


d5 


-/  iC'^X)-Q'dS-  f  {A^  ■VX  +  Y'^K)-Q’dn, 

Js^  Jo 

where  F  and  K  are  given  by  equations  (A.28)  and  (A,29),  respectively,  H  is  the  surface  curvature,  and 

Tn-il  +  ^f^)  . -• 

Setting  the  domain  integral  in  equation  (4.11)  to  zero  leads  to  the  adjoint  equations: 

(4.12)  A^-VA  +  y^K  =  0  in  fi. 

The  surface  integral  over  Soo  is  ehminated  in  the  same  way  as  that  described  in  Ref.  [28],  which  leads  to  a 
set  of  far-field  characteristic-based  boundary  conditions  for  the  adjoint  equations. 

The  surface  integral  over  Sa  has  to  be  ehminated  too.  The  contributions  from  {V'  ■  n)  and  {V  ■  s)  are 
cancelled  by  the  conditions 

(4.13)  Ta  =  (ax  +  :^A4)  p  +  -  (A  •  n)lH  +  Z(V  •  A)  +  2ax^^, 


(4.14) 


r2  =  Z^^-b(A.s)/xF-A4r^-FM  - /f(A  •  s) 


ds  '  ^  dn  '  ds  ^ 

The  terms  with  {VipE)'  ■  n)  and  (Vp'  •  n)  are  ehminated  by  the  relation 
(4.15)  Tz  =  -A4  (7^^)  . 

The  contributions  from  {pE)',  and  p'  are  set  equal  to  zero  by  satisfying  the  conditions 

&ip 

W 

dip 

dTw’ 

(4.18) 


(4.16) 

(4.17) 


An  = 

A  •  s  =  — 
( 

VA4  ♦  n  —  0. 
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These  may  be  considered  as  corresponding  to  the  no-slip  and  adiabatic  wall  boundary  conditions  (3.22).  The 
term  with  can  be  ehminated  by  the  condition 

(4.19)  A-n  =  0. 


This,  however,  conflicts  with  equation  (4.16).  This  problem  is  circumvented  by  introducing  a  term  with  Tn, 


into  ip  of  equation  (4.1),  i.e. 


..  «  .d(V  n) 

Tn-{l  +  2tl)  Q-—, 


1p  =  1p{p,Tyj,Tn,9), 


SO  that  equation  (4.8)  is  modified  to 


(4.20) 


djJ  J 


dS. 


The  associated  terms  in  equation  (4.11)  are  replaced  by  the  above  expression  appropriately,  and  equa- 
tion  (4.19)  is  replaced  by 


(4.21) 


A- 


n  =  — 


&ip 

dTn 


This  can  be  made  compatible  with  equation  (4.16)  by  imposing  the  condition 


d'tp  dip 

('■■22)  =  -W 

This  means  that  for  a  weU-posed  adjoint  problem,  there  is  a  restriction  for  the  aerodynamic  functional  :F. 
The  definition  of  J-  must  include  a  term  with  which  satisfies  equation  (4.22).  Restriction  of  the  same 
nature  was  recognized  in  Ref.  [3].  In  the  present  study,  however,  equation  (4.22)  is  proposed  as  a  general 
approach  to  ensure  the  well-posedness  of  the  adjoint  problem.  One  should  also  be  aware  that  the  combination 
of  the  continuity  equation  and  the  no-slip  boundary  conditions  dictates 


d{V-n) 

On 

implying  Tn  =  0,  so  that  introducing  a  term  with  into  T,  as  suggested  above,  does  not  modify  the 
minimization  problem  of  J-. 

The  adjoint  problem  can  now  be  summarized  as  follows.  Equation  (4.12)  in  is  to  be  solved  subject 
to  a  proper  far-field  characteristic-based  boundary  condition  on  Soo  and  the  near-field  boundary  condi¬ 
tions  (4.16)-(4.18)  on  Sa-  The  resulting  vector  of  adjoint  variables  A  is  used  for  obtaining  Ti,  T2,  and  T3 
from  equations  (4.13)-(4.15). 


4.2.  The  Gradient  Equation.  After  solving  the  flow  and  adjoint  equations,  providing  the  values  of 
Q,  A  and  T,  the  variation  of  £  becomes 


(4.23)  S£  =  SCe 

Since  0  is  a  parameter  that  describes  the  shape  of  Sa,  which  is  part  of  the  flow  domain  boundary,  the 
variation  60  implies  also  a  variation  of  the  flow  domain  fl.  As  a  result  of  this,  and  recognizing  that 


Q  =  Q(x),  X  €  n, 
A  =  A(x),  X  e  n, 
T  =  T(x),  X  €  Sa, 
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the  variation  of  fi  also  implies  a  variation  of  Q,  A,  and  X  in  the  form  of,  respectively, 


(4.24) 

/  dx 

[m 

<50), 

xe  Q, 

(4.25) 

^dx 
^d$  ■ 

sey 

xe  Q, 

(4.26) 

II 

G 

/  dx 

■  se) , 

,  xeSa 

This  leads  to  the  introduction  of  the  notion  of  the  deformation  velocity  Co  (Refs.  [8],  [13],  [28]): 

(4.27)  w(x)  =  X  •  <50, 

where 

X-iXx  Xy)  QgJ 

The  Cartesian  components  of  w  are  defined  as 

Wx  =  X®  • 

^V=Xy-  se. 

The  normal  and  tangential  components  of  w  are  written  as 

Wn  =  Xn  • 

'^s  =  Xs- 


where 

(xsJ  \%  -nxjKXyJ 

Expressions  (4.24)-(4.26)  can  now  be  written  in  the  form 


=  VQ  •  oj,  X  G  fl, 
=  VA-a5,  xeQ, 


xeSa- 


These  represent  the  so-called  convective  variations  of  Q,  A  and  Y,  respectively.  The  convective  variation 
refers  to  the  interpretation  that,  the  domain  moves  in  space  with  the  speed  w  which  gives  rise  to  the 
deformation  of  the  domain  boundaries. 

As  a  complement  to  the  convective  variation,  the  notion  local  variation  can  be  introduced  for  the  vari¬ 
ations  Q',  A',  and  X'.  The  local  and  convective  variations  constitute  the  total  variation  represented  by  a 


material  derivative: 


Q  =  Q'  -h  VQ  •  w. 


A- 

X- 


A'  -I-  VA  •  LJ, 


Ws, 
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where  the  first  and  second  terms  in  the  right-hand  sides  are  the  local  and  convective  variations,  respectively. 

The  concept  of  material  derivative  can  also  be  applied  to  the  variation  of  geometric  properties,  the 
formulae  of  which  are  given  in  Ref.  [28]  as 

•  The  material  derivative  of  a  unit  normal  vector  n; 

(4.28)  + 


where  H  is  the  surface  curvature. 

•  The  material  derivative  of  a  unit  tangential  vector  s: 


(4.29) 

The  material  derivative 

(4.30) 


of  a  surface  element  ds: 

dS  =  -I-  dS. 


•  The  material  derivative  of  a  volume  element  dCl: 


(4.31) 


dCl  =  (y  • 


For  functionals  of  the  form 


$1  = 


use  can  be  made  of  the  material  derivative  formula: 


(4.32) 


$1  =  /  fdn  -  f  fcJn  dS, 
Jn  Js 


whereas  for  functionals  of  the  form 


the  material  derivative  can  be  written  as 

dS. 

These  material  derivative  formulae  are  applied  for  obtaining  the  variation  of  SC.  Equations  (4.33)  and  (4.32) 
are  applied  for  each  surface  and  domain  integral,  respectively,  which  appears  in  the  expression  of  C.  It  is 
noted  that  the  terms  with  in  formulae  (4.33)  and  (4.32)  must  be  disregarded,  because  these  refer  to  the 
local  variations  which  have  been  eliminated  by  the  solutions  of  the  flow  and  adjoint  problems. 

The  far-field  boundary  and  the  trailing  edge  can  be  assumed  fixed  with  w  =  0,  and  the  corresponding 
terms  can  be  dropped.  For  the  sake  of  brevity,  a  tilde  notation  is  introduced  for  the  convective  variations. 
The  variation  of  C  due  to  the  variation  of  B  can  be  derived  from  equation  (4.6).  With  the  no-slip  and 


(4.33) 


$2=  / 
Js 


f  +  ff-uj  +  f(y^  +  HLj, 
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adiabatic  wall  boundary  conditions  (3.22)  taken  into  account,  this  leads  to 


(4.34) 


+  -  (A  •  n)p  +  (A  •  ^Tu,  j  (  +  HuJr 


+(A  •  n)fn  -  (A  •  n)(7  -  l)pE  -  (A  •  n  +  A  •  n)p 


+(A  •  +  (A  ■  5  +  A  •  -  |^Ai  +  p{V  •  n) 


-\-X4T^(V  •  ^  +  A47— (Ve  '  n  +  Ve  ♦  n)  +  (F  •  VA)^;^^ 


+Ti{V  •  n)  +  r2(y  •  s)  +  r3(Ve  •  n  +  Ve  •  n)J  dS. 

where  H  and  s  are  given  by  equation  (4.28)  and  (4.29),  respectively,  while 

pE  =  '^{pE)  -dj, 

'e  =  Ve-uj, 

^^(Vu.u}\ 

[Vv  ■  (3  J’ 

■=_  /VA2  -oiA 

“VVAa-wj’ 

=  (VV  ■n)-s  +  (VV  •  ^)  •  s  +  (VV  -nj-'s, 

Tn  =  (VV  ■  n)  •  n  +  (VV  •  n)  •  n  +  (VV  -n)  -n. 

Substituting  the  expressions  for  Ti,  T2,  and  T3  given  by  equations  (4.13)-(4.15),  into  equation  (4.34)  yields 

f  r^/..  ,  (^\  xa 


sc=f 

JSa  l^P 


+  n)p  +  (A  •  +  Hu); 


+(A  •  n)Tn  -  (A  •  n)(7  -  l)pE  -  (A  •  n  +  A  •  n)p 
+(A  •  s)fiu  +  (A  •  s  +  A  •  k)Tw  +  (F  •  VA)w„ 

+  Z(V.A)  +  2/.^i^V^.n)  +  /x  ^i^  +  ^i^-ff(A.s)'j(^.J)  dS. 


The  adjoint  boundary  conditions  (4.16)-(4.18)  and  the  condition  (4.22)  cancels  the  contributions  from  pE^ 
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and  Tn,  so  that  the  above  equation  reduces  to 


+  HlOn 


—(X  •  n  +  A  •  n)p  +  (A  •  s  +  A  •  +  (F  •  VA)cu„ 


+  ,(v.X)  +  2,3^)(C.»-)  +  m  ^  +  is. 


d{\-s)  a(An) 

dn  ds 


Expanding  A,  V  and  F  •  VA  gives 


(^  -  {X-  n)p  +  (A  •  j 


’  +  HuJa 


—  I  (VA2  *  (j^)  '^x  +  (VA3  •  a;)  riy 


+  HlOs]{X^^]p 


+  f  (VA2  '  (jij)ny  —  (VA3  •  0?)  Tlx  +  f  +  HcJsj  (A  •  n)  J  Tu 


8X2  ,  ^ 

aa:  +  dy 


~  (  "^xx  ~cr~  "b  (  Q_  "I  Q~"  "b  '^yy  q..  9  '  VA4 


a^ 

dy 


+  Z(V  •  A)  “t  2//- 


a(A-n) 


a(A-s-)  ,  a(A-n) 

^  dn  ds 


^  ^(Vu  ■Cj)nx  +  (Vv  ■dj)nyj 
-  H(X  ■  s)'|  ^(Vu ■Co)ny  -  (Vu  ■(3)nx^ 


The  gradient  of  the  aerodynamic  functional  f  with  respect  to  the  design  variables  0  can  be  obtained  from 


d0  ~  60' 


which  can  be  elaborated  by  using  the  definition  of  x,  equation  (4.27),  and  equations  (A.11)-(A.13)  to  give 

f  1  [  (m)  +  (’» -  o  ^  + o'  ■  + "*«) 


-( (VA2  •  x)  n*  +(VA3-x)  ny  "  (  )  P 


(VA2  •x)’^y  -  (VA3  •  x)nx  +  (  +'^A:s  j  (A -n)  j  Tu 


(dX2  aAsA  ,  ,  2  2^(9Xs^dX2\\ 

\2nxnyl^^^  QyJ  +  i%  ^  V  aa;  ay  j  J 


-(^•VA4)x„+  l(V-A)  +  2y 


a(A  •  ft) 


((Vu  ■x)nx  +  (Vu  ■x)ny) 


d(X  •  ^  0(A  •  n) 


9n  05 


-  if  (A  •  s)  (^(Vu  •  x)  -  (Vu  •  x)  ) 
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5.  Numerical  Procedure  and  Computational  Results.  The  RANS  equations  (3.3)  are  solved  by 
means  of  HI-TASK  [12],  [7]  with  the  Baldwin-Lomax  turbulence  model  [4]  implemented.  The  numerical 
procedure  in  HI-TASK  deals  with  the  time-dependent  RANS  equations  which  are  integrated  exphcitly  using 
a  five-stage  Runge-Kutta  scheme  towards  the  steady-state.  The  spatial  discretization  employs  a  cell-vertex 
finite- volume  scheme.  Jameson’s  type  of  artificial  dissipation  is  introduced  consisting  of  2-nd  order  and  4-th 
order  terms.  The  convergence  towards  the  steady-state  is  accelerated  by  means  of  a  multigrid  procedure. 
Characteristic-based  boundary  conditions  are  applied  on  Soo- 

The  adjoint  solver  employs  a  similar  numerical  scheme  as  that  used  in  the  flow  solver.  The  procedure 
deals  with  the  time-dependent  form  of  the  adjoint  equations  which  are  integrated  exphcitly  towards  the 
steady-state  using  the  same  five-stage  Runge-Kutta  scheme.  Artificial  dissipation  is  introduced  for  the 
adjoint  equations  in  the  same  way  as  that  for  the  flow  equations.  The  convergence  towards  the  steady-state 
is  accelerated  by  the  same  multigrid  procedure.  Characteristic-based  boundary  conditions  are  also  applied 
on  5oo  (Ref.  [28]). 

The  optimization  routine  FSQP  [29]  and  the  flow  solver  HI-TASK  are  integrated  with  the  adjoint  solver 
and  the  gradient  evaluator,  which  forms  the  design  code.  A  design  test  case  is  defined  representing  a 
reconstruction-type  inverse  problem.  The  target  pressure  coefficient  Cp  is  obtained  from  a  flow  analysis  of  a 
best-fit  of  the  RAE  2822  airfoil  with  the  flow  condition: 

M  =  0.73,  a  =  2°,  Re  =  6.5  x  10®, 

where  M,  a  and  Re  are  the  Mach  number,  angle  of  attack,  and  Reynolds  number,  respectively.  The  target 
Cp  distribution  is  defined  on  the  airfoil  chord  with  proper  distinction  between  the  lower  and  upper  surface 
of  the  airfoil.  The  NACA  0012  airfoil  is  used  as  the  starting  airfoil  geometry. 

The  objective  functional  to  be  minimized  has  the  form 

^  ~  2  j  ~  2  y  ~ 

where  x  is  coincident  with  the  airfoil  chord,  and  Cp^t  is  the  target  value.  The  subscripts  I  and  u  refer  to  the 
lower  and  upper  surface,  respectively.  5"  can  also  be  expressed  in  the  form 

{Cp-Cp,t)^\ny\dS. 

"  J  Sa 

The  function  ip  in  this  case  is  defined  as 

It  is  noted  that 

Cp  =  2(p-poo), 

with  p,Poo  non-diinensionahzed  hy  pooV^,  For  this  case, 

=  2{Cp  —  Cp^t)\'f^y\ 

Equation  (4.22)  requires  that  T  must  be  modified  to 

f  (Cp-  Cp,tf\ny\ dS-2  [  {Cp-  Cp,t)\ny\Tn  dS. 

^  JSa  JSa 
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The  adjoint  equation  (4.12)  is  solved  subject  to  the  adjoint  boundary  conditions  (4.16)-(4.18),  which  in  this 
case  are 

A  •  71  ” 

A  •  s  =  0, 

VA4  •  n  =  0. 

In  the  present  study,  the  deformation  velocity  has  been  formulated  as  follows, 

(5.1)  Wa  =  0  on  Sa, 

while  ujy  is  defined  by  the  curvature-continous  shape  parameterization  scheme  described  in  Ref.  [28].  This 
parameterization  scheme  has  proved  to  be  effective  in  covering  a  large  variation  of  airfoil  shape  and,  after  a 
proper  scaling,  has  shown  to  imply  an  efficient  optimization  process. 

One  purpose  of  selecting  this  test  case  is  to  investigate  the  accuracy  level  of  the  computed  gradient, 
because  the  optimal  solution  (i.e.,  the  target  airfoil)  is  known  beforehand.  The  computed  gradient  is  con¬ 
sidered  to  be  of  sufficient  accuracy  if  the  optimal  solution  can  be  obtained.  The  optimal  solution  is  assumed 
to  be  obtained  if  :F  <  10“^.  This  means  that  the  difference  between  the  actual  and  target  Cp  distributions 
is  roughly  within  0.01  (engineering  accuracy). 

Figure  5.1  shows  the  design  iteration  history.  The  engineering  accuracy  has  been  achieved  with  16  flow 
analyses.  The  optimization  was  stopped  after  the  maximum  allowable  number  of  flow  analyses  (25  analyses) 
was  exceeded.  The  Cp  distributions  and  airfoils  are  shown  in  Figure  5.2.  The  dashed  line  (the  optimization 
result)  and  the  solid  line  (the  target)  are  almost  coincident,  which  demonstrates  that  the  best-fit  of  the  RAE 
2822  has  been  closely  reconstructed. 

6.  Conclusion.  The  objective  of  the  present  study  is  to  construct  an  aerodynamic  design  methodol¬ 
ogy  using  the  variational  method  in  two-dimensional  compressible  viscous  flow  governed  by  the  Re3molds- 
Averaged  Navier-Stokes  equations.  The  focus  of  the  study  is  to  obtain  a  correct  gradient  expression. 

The  present  method  has  been  successfully  applied  for  solving  a  reconstruction-type  inverse  problem  in 
a  transonic  flow  condition.  This  means  that  the  correct  adjoint  formulation  and  gradient  expression  have 
been  obtained. 

The  numerical  result  presented  also  strongly  indicates  that  the  present  method  is  capable  of  dealing  with 
other  types  of  design  problems,  as  long  as  the  adjoint  problem  can  be  formulated  properly  as  described  in 
the  preceding  sections.  It  is  therefore  suggested  that  the  present  method  is  applied  to  more  practical  design 
cases,  such  as  those  involving  the  criteria  on  hft,  drag,  and  pitching  moment. 
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Appendix  A.  The  variations  considered  here  are  the  variations  of  the  viscous  terms  due  only  to  the  variation 
Q'  with  the  assumption  that 

•  The  variation  of  p,  with  respect  to  the  variation  of  Q  is  neglected. 

•  The  viscous  terms  on  the  far-field  boundary  are  dropped. 

The  viscous  terra  is  defined  eis 


The  variation  of  J  due  to  Q'  can  be  expressed  as 


5J=  f  \.{V-F'^)dn. 

Jci 


Integration  by  parts  yields 


(A.l) 


SJ=  -  [  A  •  (F;  •  n)  d5  -  /  P;  •  VA 
Jsa  -fn 


dVl. 
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It  is  noted  that  the  unit  normal  vector  points  toward  the  flow  domain.  Recalling  equation  (3.8)  and  intro¬ 
ducing  5Ji,  5J2,  SJz,  and  5Ji  as  Mows, 

(A. 2)  —  J  Wx  Txx  d"  '^xy)  d”  '^xy  d"  ^y  '^yy) 

(A. 3)  ^^2  ~  [  ^  '^x  {,'^xx^  d"  '^xy'^  d"  T^x'^J'  d“  TxyV  Qa:) 

d-  riy  {r'^yU  ■+  T^iyV  +  Txyu'  -f-  Tyyv'  -  q'y)  dS, 

f  f  d\o  ,  5A2  /  ^^3  f  __ 

(A.5)  5Ja  =  /  {r'xx'ti  d-  TxxU  4-  r^yV  -I-  TxyV  -  9x) 

*/  O  L 

^^4  ■ 

+  '^yy'^  +  '^yy'^  ~  Qy)  Qy 

equation  (A.l)  can  be  written  as 

(A.6)  =  SJi  —  SJ2  — 

The  variation  8Ji  will  be  dealt  with  first.  Introducing  a  local  coordinate  system  (n,s),  where  n  and  s  are 
coincident  with  the  local  normal  and  tangential  direction  on  the  surface,  respectively,  the  component  of  the 
viscous  stress  tensor  can  be  written  as 

(A.7)  Txx  =  l  +  -n^+li^nyHinyiV  .fi)-  nx{V  .^) 


^  ,  ^+r'  ^-1-/  dfi, 

dx  dy  dx  yy  dy  \ 


3^4 

xU  +  Txxu'  4-  r'xyV  4-  Txyrf  -  q'x) 


+2^  + 
d-^/i  n-x  4-  Uy 


+  2/i.  Tlx  ^y 


_  d{V.n)  d{V.^ 

~  dn  ^  ds 


+2a  + 

+  ZIJ,  Uy  +  Ux  „ 


dn 

ds  j 

'.X  (V -71)  + 

n,(f  •!)) 

d{V-^ 

d{V-n)\ 

dn 

ds  J 

9(f.s)  diV-n)'' 

"-y  On  '  dn  J  dn  ds  J 

d(V-n)  d{V-^\  ,2  2n  d{V-^  ,  d{V  n)\ 

(A.9)  Txy-^unxUy  j+n{ny  n^)  ^  4-  j 

-2/x nx  riy  H{V  ■n)-i2{nl  -  nl  )H{V  •  I). 

where  H  is  the  surface  curvature.  The  variation  can  be  worked  out  by  using  the  expression  for  Txxj  '^yyi 
and  Txy  given  above.  After  some  algebraic  manipulation,  one  obtains 


SJi=  f  (A.n)  (Z  +  2/i) 
JSa 


d{r-n) 

dn  ds 


4-(A  •  n)lH{V'  ■  n)  -  (A  •  s)^H{V'  •  s)  dS, 
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where  A  is  an  adjoint  velocity  vector  with  A2  and  A3  being  the  Cartesian  components.  The  terms  with  the 
tangential  derivative  can  be  worked  out  with  integration  by  parts  to  yield 


6J1  =  (A  •  n)l{V'  •  +  (A  •  s)tx{V'  •  n) 

+  /  [(A-n)<  +  (A-3>;-  ;x^^^-(A.n)ZF)(y'.n) 

J  Sa  / 

-  +  (A  •  {V'  •  s)  I  dS, 


where 


•w  n 


dn 

+  2/a) 


d{V'  -n) 
dn 


while  u  and  I  refers  to  the  upper  and  lower  traihng  edge,  respectively.  If  the  surface  Sa  is  assumed  smooth 
(i.e.,  a  sharp  trailing  edge  is  assumed  to  have  a  large,  but  finite,  curvature),  and  A  as  well  as  V  are  continous 
accross  the  trailing  edge,  the  first  two  terms  on  the  right-hand-side  vanish,  such  that 


(A.10) 


SJi  =  J^  (A.n)r;-b(A.s>;-  -  (A  •  n)ZffVy' •  n) 


Z^^  +  (A-5)/./fJ  {V'-^ 


dS. 


To  obtain  the  expression  for  SJ2,  use  is  made  of  the  no-slip  boundary  conditions  (3.18)-(3.19)  which  imply 

d{V-n) 

dn 

ds 

diV-n) 

ds 

It  is  noted  that  the  first  equation  above  is  identical  with  the  continuity  equation  taken  on  the  airfoil  surface 
with  u,v  =  0.  Equations  (A.7)-(A.9)  can  now  be  written  as 


(A.11) 

'^xx  —  2  nx  ny  'T'm , 

(A.12) 

"”2  nx  ny  Tjxjj 

(A.13) 

'^xy  ”  (  ^y  ^x 

Substituting  these  and  equations  (3.12)-(3.13)  into  (A. 3)  gives 


Noting  that 
(A.14) 
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or  in  terms  of  the  total  energy, 
(A.15) 

the  variation  6J2  can  be  written  as 


•^d  =  d? 


{pEY  p'  \ 

P  (7-  l)pj  ’ 


f  Xi  TUV'-d)  +  a^ 

Js. 


(A.I6)  6J2  =  A4  j  ] 

Next,  the  variation  of  is  obtained  as  follows.  Equations  (3.9)-(3.10)  are  substituted  into  equation  (A.4) 
which  yields 


( dv!  5A2  dv^  dAsA  / 5A2  dX3\  ( dv^  du' 

^  ^\dx  dx  ^  dy  dy  )  ^\dy  0a;  /  \ dx  dy 


Rearranging  the  coefficients  of  (i  gives 


8X2  du'  8X2  du' 
8x  dx  dy  dy 


8X3  du' 
\  dx  dx  dx  dy 


8X3  dv'  8X3  dv' 
dx  dx  ^  dy  dy 

'8X2  dv'  8X3  dv' 
dx  dx  dy  dy 


This  can  be  written  in  a  compact  form  as 


5Jz  = 


=  J  {V  ■X)l{V  ■V')  +  n  VA2  •  W  +  VAg  •  Vu' +  ^  •  Vti' +  •  Vu' j  dn. 


Using  the  following  vector  identities, 


(V  •  A)Z(V  •  U')  =  V  •  (i(V  •  A)U')  -  V(i(V  ■  A))  •  V', 
fJ.VX2  •  Vtt'  =  V  •  /iu'VAg  -  u'(V  ■  /iVAg), 

/iVAg  •  W  =  V  •  /^u'VAg  -  u'(V  •  /iVAg), 

dX  ^  ^  ,8X  ,  -  dX\ 


dx  -  ,  ^  ,8X  ,  -  dx\ 

'^dy  dy  ^dy ) 


one  obtains 


(A.  17)  5J3  =  j  \^-{l{V-X)V')  +  V-yu''VX2  +  W-yv'VX3  +  V-pu'^  +  V-nv'^ 

-L\ 


+  ^i^^  +  v-MVAg  +  v-Ml^U'  <m. 


20 


The  terms  having  the  divergence  form,  which  are  collected  under  the  first  integral,  are  worked  out  using  the 
Gauss  theorem  (with  the  integral  over  Soo  neglected): 


[  V  •  (i(v  •  X)V')  dVt  = 

Jo. 

I  V  •  ^u'V\2  dVt  — 
Jet 

I  V  •  /iu'VAs  dQ  = 
Jet 

fdX 

I  V  •  au'  77”  dO.  = 
Ja  dx 

[  v-fj,v'^dn  = 

Jn  9y 


-f  1{V  .X){V' ■n)dS, 
JSa 

—  /  iJbu^V\2  ’^dS^ 

JSa 

—  I  juu' VA3  •  n  dS, 

JSa 

—  /  au  —  -ndSj 

JSa 

f 

Js/''  dy 


•  fi  dS, 


Now,  the  following  notations  are  introduced 
(A.18) 

(A.  19)  r 

(A.20) 


-  -  8X2 

^xx  =  ZV  •  A  + 


yy 


lV-X  +  2y,^, 
dy 


r  - 

^  dyj- 


which  may  be  considered  as  the  elements  of  an  adjoint  stress  tensor  because  of  the  close  resemblance  with 
Txxf  T'xyj  and  Tyy.  After  some  manipulations,  equation  (A.  17)  can  be  written  as 


(A.21) 


SJ3 


=  ^  f  /(V  .  A)  + 
JSa 


2fJL 


d{X  •  n) 


+/X 


dn  y 
0(A  •  ^  9(A  •  ft) 


dn 


+ 


ds 


(V'-ft) 


H{X-^ 


dS 


f  r  / 

Jn  I  \  9x 


ar* 


XX  xy  I  ^1 


dy 


dVa^y  .  dT. 


dx 


dy 


yv  >  y> 


dJ^. 


The  variation  6Ji  (A.5)  is  worked  out  using  equations  (3.9)-(3.17).  This  gives 

„  /  au'aA4  au'aA4\  .  fdv'  ,  dv/ 

^  dx  dx  dy  dy  )  ^\dx  dy  /  \ 

f  dX^  dXiX  .  (  dXi  ^  dXA  , 

^  ay  J  “  +  dx  ^  ay  j  ^ 

\  aa;  ax  dydy) 


dXi  dXi 


+7 


J± 

Pr 
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After  some  manipulations,  one  obtains 


SJa  =  l\{y-  VA4)1(V  •  y') 
Ja  . 


/  du'  d\i  du'  8X4 \ 

dx  dx  ^  dy  dy  J 


- - -  - -  ,  ,  81/8X4  .  dv'dX4 

n..  I  dx  8x  8y  dy 


f  du' 

du''\ 

I  aA4 

/  dv' 

dv'X 

van 

dy  J 

van 

+  7-5 


9A4 
'  dx 


8X4 


dy 


■f  X'xy  Q  1  ^  "b  (  Txy  b  "rj 


aA4 


y  f  8e'  8X4  de'  8X4 


\dx  dx  dy  dy 


dx 

dCt. 


‘yy 


8X4 

dy 


A  compact  form  can  be  written  as  follows, 


6J4 


Jn 


{V  ■  VA4)/(V  ■  V') 

+yuVX4  ■  Vn'  +  yvVX4  ■  Vv'  +  y^V  •  Vu'  +  y^V  ■  Vv' 
(  8X4  dX4\  ,  ,  f  8X4  ,  8X4\  , 

+  j  ^  ^  dy)'’ 


+7^^  (ve'.VA4 


dfl. 


The  following  vector  identities  are  considered, 


{V  .  VA4)KV  •  F')  =  V  •  {1{V  .  VA4)f ')  -  V' .  Vl{V .  VA4), 
fjbuVX4  *  Vu'  =  V  •  /itfcw  V  A4  t//(V  •  /inVA4), 

//VVA4  •  Vi;'  =  V  •  iJ>vv'WX4  —  i;'(V  •  pivVX^)^ 


7^(Ve'  •  VA4)  -  :^(V  •  Me'VA4  -  e'(V  •  MVA4)), 
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for  obtaining 


8Ja=  f  V  ■  {1{V  ■WX4)V')  +  V  ■  fiuu'VXi  +  V  ■  fivv'VXi 

Ja  L 

+V-/i^«'y  +  V-/i^u'F+^V-Me'VA4  dn 

^  dx  dy  Pr 

L  [  (’■"^  It  -  ^  ^  ^  “' 

+  -  |-(,r .  VA,)  -  V  .  ..VA,  -  V  .  v)  «' 


--f(V-^VA4)e'  do. 

Pr  J 

Applying  the  Gauss  theorem  and  the  no-slip  boundary  condition  for  the  first  integral,  and  introducing 
and  as 


(A.22) 


(A.23) 


(A.24) 


^  ,  9A4  ,  9A4 

=  {I  +  ’ 

/,  ^  .  9X4  ,  8X4 

^yy  —  {I  +  2fj)v  +lu  , 

f  8X4  8X4\ 

Qy  dx)^ 


(A.25) 


SJ4  =  -f  7^(VA4-n)e'dS' 
JSa 


f  f  d\^  5A4  XX  d^xy  ^  f 


5A4  9A4  d^xy  9^yy  A  / 

dy  dx  dy  J 


-i(V.MVA4)e'J  da. 

Substituting  equation  (A.15)  into  e'  in  the  surface  integral  and  equation  (A.14)  into  e'  in  the  domain  integral 
results  in 


(A.26) 


P’ 

(7-  l)p 


dA4  dX4  d^xx  d^xy\  , 

dx  dy  dx  dy  ) 


dX4  .  aA4  d'^xy  V  , 


I  '^yy 


dx  dy  dx  dy 


a?{V-yiWX4)\  fp'  f/ 


(7-l)Pr  )  \p  p 
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Substitution  of  equations  (A.IO),  (A.16),  (A.21),  and  (A.26)  into  (A.6)  yields 

(A.27)  |^(A  •  +  (X  •  n)r; 


-  (A  •  n)lH  ]  (V'  .n)  -  +  (A  •  |  (V'  •  s) 


,^(A-n) 

ds 


■^XA'^wiy'  ■  ^  "i"  A4tt 


2 /X  W{pE)' -n  ypf  ■  n 
'  Pr  p  (7  -  l)p 


/  SPxx  dVxy  ^A4  ^A4  d^xx  .  d"^xy\  j 

_  V  dy  dx  dy  dx  dy  )  ' 

,  (dVxy  ,  dVyy  _  aA4  ,  aA4  ,  5^.^  , 

+  +  "^’^92/+  aa;  +  dy 


dx  dy  dx  dy 


a2(V-/xVA4)\  /p'  p' 


+  -  y  dn. 

(7-l)Pr  J  \p  p) 

The  domain  integral  can  be  expressed  in  terms  of  the  conservative  flow  variables  by  using  the  transformation 

U'  =  FQ', 

where  F  is  the  Jacobian  of  the  primitive  flow  variables,  U  =  (p  u  v  py,  with  respect  to  Q, 


(A.28) 


—  ±00 
P  P 

-  0  J  ° 

{'y  -  l)(u^  +  ^  I  ^ 

— - -  -{y-l)u  -(7-l)u  7-1 


The  coefficients  of  U'  in  equation  (A.27)  can  be  collected  into  a  vector  K  defined  as 

/  a2(V  •  PVA4)  1  \ 


(7  -  l)Pr  p 


(A.29) 


OYxx  . 

dFxy 

1 

1  ^ 

1 

dXi 

-1- 

d'^xx 

4- 

a^*y 

dx 

dy 

•XX  Q 

OX 

dy 

1 

dx 

dy 

dTxy 

dFyy 

dXi 

dXi 

4- 

d^xy 

_J_ 

d%y 

dx 

dy 

dx 

dy 

dx 

dy 

a^(V  •  PVA4)  1 
(7-l)Pr  p 
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Finally,  equation  (A.27)  is  written  as 


(A.30)  [(^  •  ^yn  + 


-  fJ' 


^  -  (A  •  nw]  {V'  •  n)  - 


ds 


J 


ds 


..  ::a  ,  A  _2M  'VipEY-n  Vp'  n 

+A4T^(F  •  s)  +  A4a  —  - 


-  Z(V-A)  +  2/x 


0(A  •  n) 
dn 


(F'-n)-/. 


0(A  •  s)  a(A  •  n) 


9n 


5s 


-a^^(VA4-n)^ 


( jpE)'  p' 

P  (7  - 1)/> 


dS 


+  f  Y^K  Q'dn. 

Ja 


iv'-^ 


F(A.s)j  (F'-s) 
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number  of  aerodynamic  analyses 


Fig.  5.1.  Optimization  histoTy.  JVf  =  0.73,  a  —  2°,  Rc  =  6.5  x  10^. 


Fig.  5.2.  Cp  distribution  and  airfoil  geometry.  M  =  0.73,  a  =  2°,  Rc  =  6.5  x  10*^ 


